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Abstract 



Our understanding of the mechanism by which topological defects are 
formed in symmetry breaking phase transitions has recently changed. We 
examine the non-equilibrium dynamics of defect formation for weakly-coupled 
global 0(N) theories possessing vortices (strings) and monopoles. It is seen 
that, as domains form and grow, defects are swept along on their boundaries 
at a density of about one defect per coherence area (strings) or per coherence 
volume (monopoles). 

I. INTRODUCTION 

The formation of topological defects during symmetry breaking phase transitions is 
generic to many physical systems. In particular we cite the vortices and monopoles of 
superfluid ^He and ^He and the vortices (flux-tubes) of high- and low-Tc superconductors. 
Similar defects, cosmic strings or monopoles, most likely appeared in the early universe at 
the GUT-scale phase transition. All of these systems are described by some form of quantum 
field theory and, due to the phase transition, their dynamics is intrinsically non-equilibrium. 
They therefore provide a good means to test non-equilibrium field-theory experimentally 
over a wide range of energies. 

Roughly, the dynamics of defect formation proceeds as follows From some initial 
state, which is not too far from thermal equilibrium, some change in the bulk properties of 
the system, such as pressure or volume, induces a phase transition. During this transition, 
the scalar fields which describe the order parameter fall from the false vacuum into the 
true vacuum, choosing a point on the vacuum manifold at each point in space, subject to 
the constraint that they must be continuous and single- valued. We shall limit ourselves to 
weakly first order or continuous transitions, for which this collapse to the true vacuum occurs 
by spinodal decomposition or phase separation. The resulting field configuration is one of 
domains within each of which the scalar fields have relaxed to a constant vacuum value. If 
the theory permits defects, it will sometimes happen that the requirements of continuity 
and single valuedness force the fields to remain in the false vacuum between some of the 
domains. For example, in the case of a complex scalar field producing vortices, the phase 
of the field may change by an integer multiple of 2n on going round a loop in space. This 
requires at least one zero of the field within the loop, which signifies the presence of a region 
of unbroken phase. Each zero has topological stability and characterises a vortex passing 
through the loop. When the phase transition is complete and there is no longer sufficient 
thermal energy available for the field to fiuctuate into the false vacuum, the topological 
defects are frozen into the field. From then on, the defect density alters almost entirely by 
interactions of defects amongst themselves, rather than by fiuctuations in the fields, see for 
example 0. 

The major question then, is what fixes the initial defect density and the defect correla- 
tions? Only then can the subsequent evolution of defect networks be determined with any 
accuracy. It was first argued that topological defects should be frozen in at the Ginzberg 
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temperature Tq [Q], the temperature above which there is sufficient thermal energy available 
for the field to fiuctuate into the false vacuum without cost . If so, the defect number 

would be strongly fiuctuating above the Ginzberg temperature but frozen in below it. In 
this case, the relevant scale for the initial defect density would be the coherence length, 
^{Tg), of the Higgs field (or fields) at the Ginzberg temperature. For example, in vortex 
production for the U{1) theory mentioned above, the initial vortex density (i.e. the number 
of vortices passing through unit area) would be k/^'^{Tg), where k is a constant of order 
unity. Similarly, in monopole production the monopole density would be expected to be 
k/^^{Tg), for similar n. Thereafter, defect forces are assumed to take over. 

Recently, however, more compelling pictures of the way in which the initial density of 
topological defects is fixed have been proposed. While the mechanism outlined initially is 
almost certainly correct, in general it is unlikely that the Ginzberg temperature is relevant 
to anything other than a thermally produced population of defects. For the cases of interest, 
for example an expanding universe, we expect that, as the system is driven from some initial 
thermal state towards the phase transition there comes a point when the rate at which the 
transition is driven is too fast for the evolution of the field to keep up The transition 
may now be viewed as a quench and it is not clear that either temperature or free energy 
mean anything at all. At this point any defects within the field are assumed to be frozen 
in until the transition is complete. Upon completion, the field will try to return to thermal 
equilibrium. At sufficiently low temperatures, however, the return to equilibrium by thermal 
processes will be so slow that the evolution of the initial defect density thus produced will 
be almost entirely by interactions between the defects themselves. 

Thus, in this scenario, the vortices cease to be produced, not at the Ginzberg tempera- 
ture, but when the scalar fields go out of equilibrium The relevant scale which determines 
the defect density is the coherence length, ^(t) at this time, and for some time onwards, 
rather than the coherence length at the Ginzberg temperature. The only remaining uncer- 
tainty is how good an approxmiation it is to say that the defect number is frozen into the 
field from the time when it first goes out of equilibrium until the time when it lies in the 
vacuum manifold almost everywhere and its evolution can be viewed as the interaction of 
the defects. It is this question which we address here. 

II. THE MODEL 

In the following we consider a class of theories where the broken and unbroken symmetries 
are global, thereby guaranteeing that they pass through a second order transition. [Had they 
passed through a strongly first order transition, the mechanism for the transition, bubble 
nucleation, would lead to different consequences from those outlined below, although it 
might still be a good approximation to say that the defect density is frozen in when the 
field first goes out of thermal equilibrium]. We assume that the change of symmetry is 
sufficiently rapid that the fields are unable to respond immediately, but evolve by means of 



Although thermal equilibrium is mode dependent, this does not matter for the crude argument 
repeated here. 
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phase separation or spinodal decomposition and domain formation. 

We shall consider the simplest theory, one of massive relativistic scalar fields 0^, where 
a = 1, . . . , N, in D spatial dimensions, transforming as the fundamental representation of 
a globally 0{N) invariant theory. Changes in the environment cause the symmetry to be 
broken to 0{N — 1) (i.e. as given by the generalised 'wine-bottle' potential) leading to a 
theory of one massive Higgs boson and — 1 massless Goldstone bosons, with the vacuum 
manifold S^~^. Since the nth homotopy group n„ of the n-sphere is n„(5'") = Z, the group 
of integers, the theory possesses global monopoles if N = D and global strings if N = D — 1. 
We are primarily interested in D = 3 dimensions, for which the 0(3) theory possesses 
monopoles, and the 0(2) theory possesses strings. [However, the vortex production in the 
D = 2 Kosterlitz-Thouless transition has some interest, al though we shall not pursue it 
here] . 

The transition is realised by the changing environment inducing an explicit time- 
dependence in the field parameters. Although we have the early universe in mind, we 
remain as simple as possible, in fiat space-time with the 0-field action:- 

The t-dependence of m?{t) and \{t) is assumed given and adjacent 0{N) indices are summed 
over. 

We wish to calculate the evolution of the defect density during the fall from the false 
vacuum to the true vacuum after a rapid quench from an initial thermal state. The simplest 
assumption, which we shall adopt, is that the symmetry breaking occurs at time t = to, 
with the sign of m?{t) changing from positive to negative at to- Further, after some short 
period At = t — to > 0, m?{t) and A(t) have relaxed to their final values, denoted by m? 
and A respectively. The field begins to respond to the symmetry-breaking at t = to but we 
assume that its response time is greater than At, again ignoring any mode dependence. 

To follow the evolution of the defect density during the fall off the hill involves two 
problems. The first is how to count the defects and the second is how to follow the evolution 
of the quantum field. We take these in turn. 

III. COUNTING THE DEFECT DENSITY 

To calculate the defect density requires knowledge of Pt[^], the probability that, 
the measurement of the field 0(t, x) = (0i, 02, . . . , ^at) would yield the result $(x) = 
($1, $2, • • • 5 '^'Af)- This is obviously a consequence of both the initial conditions and the 
subsequent dynamics. We follow Halperin |]^ in adopting a Gaussian distribution for the 
field of the form:- 

=Arexp(^--l|d°xd°y $„(x)i^„6(x-y;t)$6(y)), 

with Kab = SabK and JV a normalisation. The circumstances under which a Gaussian is 
valid will be examined later. For weakly coupled theories we shall see that, for short times 
after to at least, a Gaussian pt[$] will occur. If this is taken for granted it is relatively 
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straightforward to calculate the number density of defects. Postponing the calculation of K 
until then, we quote those of Halperin's results that are relevant to us. 

Suppose that the field 0(t,x) takes the particular value $(x). We count the vortices 
by identifying them with its zeroes. The only way for a zero to occur with significant 
probability is at the centre of a topological defect so, but for a set of measure zero, all zeroes 
are topological defects. 

Consider the 0{D) theory in D spatial dimensions, with global monopoles. Although less 
relevant than strings for the early universe they are slightly easier to perform calculations 
for. Almost everywhere, monopoles occur at the zeroes of $(x), labelled Xi,i = 1,2, . . ., at 
which the orientation $(x)/|$(x)| is ill-defined. A topological winding number rii = ±1 can 
be associated with each zero Xi by the rule:- 

rii = sgndet(9a$fc 

Monopoles with higher winding number are understood as multiple zeroes of $(x) at which 
the Hi are additive. The net monopole density is then given by:- 

Pnet (x) = ^ ni 5(X - Xi) . 

i 

The volume integral of this gives the number of monopoles minus the number of anti- 
monopoles. The correlations of pnet give us information on monopole- (anti)monopole corre- 
lations but, in the first instance, we are interested in the cruder grand totals. The quantity 
of greater relevance to us, is the total monopole density:- 

P(x) = ^5(X- Xi), 
i 

whose volume integral gives the total number of monopoles plus antimonopoles. in the 
volume of integration. 

Now consider an ensemble of systems in which the fields $ are distributed according to 
at time t. Then, on average, the total monopole density is:- 



(p(t)) = /E5(x-Xi 

i 

where the triangular brackets denote averaging with respect to Pt[$]. That is:- 

(F[<l>])i = / P$F[<l>]pt[<l>], 



with pt[^] normalised so that / T'$pt[$] = 1. The translational invariance of the Gaussian 
kernel of the probabililty density ensures that p(t) is translationally invariant. 
In terms of the fields vanishing at Xi, p(t) can be re-expressed as:- 



{p{t)) = (5^[$c(x)] \det{dM^))\),. 



^This counting procedure alone gives no information about the length distribution of the defects. 
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The second term in the brackets is just the Jacobian of the transformation from x to $(x). 

It follows from the Gaussian form of the probability density that the individually and 
independently Gaussian distributed with zero mean, as 

($„(x)$6(x))t = 5„;,W^(|x-y|;t), 

where iy(|x — y|;t) = i^~^(x — y;t). So also are the first derivatives of the field da^b, which 
are independent of the field:- 

($e(x)a„$6(y))t = 

due to the fact that W is dependent only on the magnitude of x — y. 

Thus, the total defect density may be separated into two independent parts:- 

(p(t)) = (5^[$(x)]),(|det(a.$,)|),. 

The first factor is easy to calculate, the second less so. 
Consider first the delta-distribution factor:- 

{6m^o)])i = J V^5m^o)] exp|-i|(i^X(/^y$(x)ir(x-y;t)$(y)| 

^ jiajv^ exp|ia$(xo) - ^ / d^xd^y$(x)X(x - y; t)^{y) |, 

where 0{N) indices and integrals over spatial variables have been suppressed and = 
da/27r. On defining a — 5(x — xo)q;, we find:- 

(5[$(xo)])t ^ l^ajv^expy rf^x (^ia(x)$(x) -^J ci^xd^y$(x)K(x - y; i)$(y)) |<5(> 
= J ^aexp|-^y ci^xci^ya(x)W^(|x-y|;i)a(y)| 
= J (laexpj — - J (/■^X(i^y(5(x — Xo)Q;iy(|x — x'l; t)Q;(5(x' — Xq) | 



Y. 
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2^^\^/K^J 27r($$)' 27vW{0;t) 

Consider now the second factor. Writing out the determinant explicitly for N — D — 2, 
and exploiting the fact that the field is Gaussian, we have:- 



( I det{daU^)) I )l 



det((9a0b(x)) 



(9i0i9202) + (9l02920l) - 29i02920l9i0i9202 



The first term may be factorised into a product of two Gaussian variables and calculated as 
follows:- 

{{dicl>,d2h)')t = {{dM')t{{d2<l>2?)t 

= [-5ndidiW{x;t)] [-^^2^2^^(1x1;^)] = [d.diW t)]'', 
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where l^(|x|;i) is as before. Fourier transforming the two-point function, we find:- 



(0282! 



{ (9101^202)' )t = {^,8^ J e'^-^iy (k; t) 

cos\e)W{k;t) fk 

A similar result applies for the second term whereas the third term vanishes. The final result 
is:- 



(Pit)) = C^iv 



W"{0; t) 



W(0;t) 



N/2 



where the second derivative in the numerator is with respect to a; = |x|. Cat is I/tt^ for 
N = D = 3 and l/27r for N = D = 2 (had we performed the calculation in D = 2 
dimensions), the difference coming entirely from the determinant factor. 

Let us now consider the case of global strings in D = 3 dimensions that arise from the 
0(2) theory. Strings are identified with lines of zeroes of 0(t,x) = $(x) and the net vortex 
density (vortices minus antivortices) in a plane perpendicular to the i-direction is:- 

Pnet,i(x) = (5^[$(x)]eij-fc(9j$i)(9fc$2), 

with obvious generalisations to = D — 1, for all in terms of the Levi-Cevita symbol 
As before, the total vortex density is of more immediate use. On a surface 
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perpendicular to the i-direction this is:- 

Pnet,i(x) = 5'[$(X)] \eijk{dj<^i){8k<^2)\. 

in analogy with the monopole case. The expectation value of this total density, when 
calculated as before reproduces the same expression:- 



(pit)) = C 



N 



W"{0;t) 



W{0;t) 



N/2 



but for A^ = D — 1. For the case of interest, N = 2 and C2 = l/27r. Thus, whether we are 
concerned about global monopole or global string density, once we have calculated K{'x; t) 
we can find the total string density. Similar results apply for the correlations between net 
densities, which are important in determining the subsequent evolution of the defect network. 



IV. EVOLUTION OF THE QUANTUM FIELD 



During a symmetry breaking phase transition, the dynamics of the quantum field is 
intrinsically non-equilibrium. The normal techniques of equilibrium thermal field theory 
are therefore inapplicable. Out of equilibrium, one typically proceeds using a functional 
Schrodinger equation or using the closed time path formalism of Mahanthappa, Schwinger 
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and Keldysh [p|-p!0[] . Here, we employ the latter, following closely the work of Boyanovsky,de 



Vega and coauthors [|ll| , [r2 



Take t = to as our starting time. Suppose that, at this time, the system is in a pure 
state, in which the measurement of would give $o(x). That is:- 

0(to,x)|$o,to) = '^'ol'^'cto)- 

The probability pt^.[^f] that, at time tj > to, the measurement of will give the value 
is Ptfi^f] = |c/oP, where:- 



9^(*/)=*/ 
</>(to)=*o 



V(j) exp< iS[(f)] 



in which V(j) = Yla=i'^4'a and spatial labels have been suppressed. It follows that pt^l^f] 
can be written in the closed time-path form:- 



0±(to)=*o 



Instead of separately integrating 0± along the time paths to < t < tf, the integral can be 
interpreted as time-ordering of a field along the closed path C+ © C_ where = 0+ on 
C-L and 6 = 6^ on C_. 
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FIG. 1. The closed time path contour C+ © C_. 



It is convenient to extend the contour from tf to t — oo. Either or 0_ is an equally 
good candidate for the physical field, but we choose 0+:- 

t 



to C . 



FIG. 2. Extending the integration contour. 
With this choice and suitable normalisation, becomes:- 

PtA^f]= I , ^ 2)</>+P0-<5[0+(t)-$/] exp|^(5[0+]-5[0_])j, 

J(/.±(to)=*o IV /J 

where — $/] is a delta functional, imposing the constraint x) = ^/(x) for each 

X. 

The choice of a pure state at time to is too simple to be of any use. The one fixed 
condition is that we begin in a symmetric state with (0) = at time t = to. Otherwise, our 
ignorance is parametrised in the probability distribution that at time to, 4>{to, x) = $(x). If 
we allow for an initial probability distribution Pto[^] then is generahsed to:- 

PtA^f] = I P$P,J$]£^^^^^^^^P0+P0_5[0+(i) exp|i(^[0+] 

At this stage, we have to begin to make approximations. So that ptj[^f] shall be Gaussian, 
it is necessary to take -Pto[$] to be Gaussian also, with zero mean. All the cases that we 
might wish to consider are encompassed in the assumption that $ is Boltzmann distributed 
at time to at an effective temperature of To — (5q^ according to a quadratic Hamiltonian 
Hom. That is:- 

PtM = ($,to|e-''^1$,to) = / P03exp|^5'o[03]j, 
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for a corresponding action 5*0 [03], in which 03 is taken to be periodic in imaginary time with 
period /?o- We take 5*0 [03] to be quadratic in the 0{N) vector 03 as:- 

We stress that mo and /3o parametrise our uncertainty in the initial conditions. The choice 
oo corresponds to choosing the pt[^] to be determined by the ground state functional 
of Hq, for example. Whatever, the effect is to give an action 5*3 [0] in which we are in thermal 
equilibrium for t < to during which period the mass m{t) takes the constant value mo and, 
by virtue of choosing a Gaussian initial distribution, X{t) — ior t < to. 
We now have the exphcit form for pt^ [$ /] :- 

p. [$J= [v^ [ P03e^^°['^3] / r'0+$me'(^[^+]-^['^-])5[0+(i/^) 

- j^V<P,V<P+V<P_ exp|^5o[03] + ^(^[0+] - ^[0-])} 5[0+(t/) - $/], 

where the boundary condition B is (f)±{to) — (psito) — (psito — iPo)- This can be written as 
the time ordering of a single field:- 

pt^i^f] = I D0e'^-[^U[0+(t/)-$/], 

along the contour C = C+ © C_ ® C3, extended to include a third imaginary leg, where 
takes the values 0+, 0_ and 03 on C+, C_ and C3 respectively, for which Sc is 5'[0+], <S'[0_] 
and 5*0 [03]. 



1 





t 

> • 






C3 




c. 



FIG. 3. A third imaginary leg 

We stress again that although 5'o[0] may look like the quadratic part of 5'[0], its role 
is solely to encode the initial distribution of configurations $ and need have nothing to do 
with the physical action. Henceforth we drop the suffix / on and take the origin in time 
from which the evolution begins as to = 0. 

We perform one final manoeuvre with pt[^] before resorting to further approximation. 
This will enable us to avoid an ill-defined inversion of a two-point function later on. Consider 
the generating functional:- 

Z\j+,3-,h] = j^'^'P exp|i5c[0] +i j j0|. 



5o[03] = / d°+^x 
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where / jcp is a short notation for:- 



omitting spatial arguments. Then introducing aa(x) where a = 1, . . . , A^, we find:- 

= J T^oi j^Vcj) exp|zSc;[0]| exp|« j d^xa,(x)[0+(t/, x) - $(x)],| 
= jva exp|-i j a„$„|z[a,0,0], 

where a is the source a(t, x) = a(x)(5(t — t/). As with Vcp, Va denotes Hi^ ^^^a- 

We have seen that analytic progress can only be made insofar as pt{^] is itself Gaussian, 
requiring in turn that Z\a, 0, 0] be Gaussian in the source a. In order to treat the fall from 
the false into the true vacuum, at best this means adopting a self-consistent or variational 
approach i.e. a Hartree approximation or a large N-expansion. However, if we limit 
ourselves to small times t then pt[^] is genuinely Gaussian since the field has not yet felt 
the upturn of the potential. That is, we may treat the potential as an inverted parabola 
until the field begins to probe beyond the spinodal point. The length of time for which it is 
a good approximation to ignore the upturn of the potential is greatest for weakly coupled 
theories which, for the sake of calculation we assume, but physically, we expect that if the 
defect counting approximation is going to fail, then it will do so in the early part of the fall 
down the hill. 



V. EVOLUTION OF THE DEFECT DENSITY 

The onset of the phase transition at time t = is characterised by the instabilities of long 
wavelength fluctuations permitting the growth of correlations. Although the initial value of 
(0) over any volume is zero, the resulting phase separation or spinodal decomposition will 
lead to domains of constant {(p) whose boundaries will evolve so that ultimately, the average 
value of in some finite volume, will be non-zero. That is, the relativistic system has a 
non-conserved order parameter. In this sense, the model considered here is similar to those 
describing the A transition in liquid helium or transitions in a superconductors. 

Consider small amplitude fluctuations of 0a, at the top of the parabolic potential hill 
described by V^(0) = |m^(t)0^. At t < 0, rn?{t) > and, for t > 0, m^{t) < 0. However, by 
t ^ At, im?{t) and A have achieved their final values, namely — /i^ and A. Long wavelength 
fluctuations, for which |kp < — m^(t), begin to grow exponentially. If their growth rate 
Ffc —im?{t) — |kp is much slower than the rate of change of the environment which is 
causing the quench, then those long wavelength modes are unable to track the quench. For 
the case in point, this requires mAt -C 1. We take this to be the case. To exemplify the 
growth of domains and the attendant dispersal of defects, it is sufficient to take the idealised 



^In the latter case, given the relationship between and spatial dimension D = N, this corre- 
sponds to a large-dimension expansion 
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case, At = in which the change of parameters at i = is instantaneous. That is, m^(i) 
satisfies:- 



ml>0 iit<0, 
< if i > 



where for t < 0, the field is in thermal equilibrium at inverse temperature /?o- As for X{t), 
for i < 0, it has already been set to zero, so that Pto[^] be Gaussian. For small t, when 
the amplitude of the field fluctuations is small, the field has yet to experience the upturn 
of the potential and we can take A(t) = then as well. At best, this can be valid until the 
exponential growth |0| ^ /ie'** in the amplitude reaches the point of inflection |0| fti iJ,/\/X, 
that is /it ~ 0(ln(l/A)). The smaller the coupling then, the longer this approximation is 
valid. As noted earlier, it should be possible to perform more sophisticated calculations 
with the aim of evolving the defect density right through the transition. For our present 
purposes, however, the small time or Gaussian approximation is adequate. 

We are now in a position to evaluate pt[^], identify K and calculate the defect den- 
sity accordingly. Sc[<p] becomes 6'o[03] on segment C3 so setting the boundary condition 
0+(O, x) = 03(0, x) = 03(— i/3o,x) and we have:- 

S[<f>^] = J d^+^X [^{dM{d^<l>a) + ^/XV^], 

on C+. The Gaussian integrals can now be performed to give:- 

= Jvaexpl^-i yd°xaA|exp|^y^d^xd°ya„(x)G'(x-y;t,t)a(.(y)|, 

where G(x— y; t, t) is the equal time correlation, or Wightman, function with thermal bound- 
ary conditions. Because of the time evolution there is no time translation invariance in the 
double time label. As this is not simply invertible, we leave the a integration unperformed. 
The form is then a mnemonic reminding us that K^^ = G. 

In fact, there is no need to integrate the as since from the previous equation it follows 
that the characteristic functional (exp{i / Ja^a})t is directly calculable as:- 



|exp|i J Ja^a}^^ = J V^Ptm expji J 

- expji J d^xdV Ja(x)G(x - y; t)j„(y)}. 
Thus for example, the first factor in the monopole density p{t) is:- 

(5^[$(x)]), = ^|d7-exp(i$„(x)i„)^ 

= J dj exp|^j2G(0;i,i)} = [-iG{0;t,t)]-^/', 

with suitable normalisation, without having to invert G{0;t). Thus, on identifying 
—iG{-K;t,t) with l^(x, t) as defined earlier, p{t) becomes:- 
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(pit)) = Cm 



-iG{0;t,t) 



N/2 



where — zG(x; t,t) has to be calculated from the equations of motion, subject to the initial 
condtition. 

Details are given by Boyanovsky et al., and we quote their results, which give 
— 'iG(x;t, t) as the real, positive quantity:- 



-iG(x;t,t) 



tk.x 



coth(/5otu<(A;)/2)< 



+ 



1 + ak{cos{2u~^{k)t) — 1 



1 + Ak{cosh{2W{k)t) - 1) 



with:- 



ulik) - 
colik) 



\k\^ + mf, 



M -P 



A, 



Oik 



2^ W\k) J 



1 - 



The first term is the contribution of the unstable long wavelength modes, which relax most 
quickly; the second is that of the short wavelength stable modes which provide the noise. 
The first term will dominate for large times and even though the approximation is only 
valid for small times, there is a regime, for small couplings, in which t is large enough for 
cosh{2fit) ~ |exp(2/it) and yet /it is still smaller than the time 0(lnl/A) at which the 
fluctuations sample the deviation from a parabolic hill. In these circumstances the integral 
at time t is dominated by a peak in the integrand k^-^^'^^i^)^ at k around kc-, where:- 



{D-r 



The effect of changing [3q is only visible in the 0(1 /fit) term. In the region |x| < yt/fi the 
integral is dominated by the saddle-point at kc, to give:- 



iG(x; t, t) = W{x; t) ^ W{0; t) exp I 



-fxx 



smc 



for D = 3, where :- 



W{0;t) ^ C 



for some C, which we don't need to know. The exponential growth of G{0] t) in t reflects the 
way the field amplitudes fall off the hill ($) = 0. It is sufficient for our purposes to retain 
D = 3 only. 
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After symmetry breaking to 0{N — 1) the mass of the Higgs is mu = \fl\x with cold 
correlation length ^(0) = ■ On identifying e~^^ as 1^ we interpret:- 



m = i8V2Y/'^ tm, 

as the size of Higgs field domains. This t^/^ growth behaviour at early times is characteristic 
of relativistic systems (with a double time derivative) with a non-conserved order parameter. 

To calculate the number density of defects at early times we have to insert this expression 
for —iG or W into the equations derived earlier. Expanding W{x; t) as:- 

W(x;t)^W{0;t)e^p(g^)(l-^^ + o( 



e{t)j V sew \m 



and substituting in () we find:- 





, - , /"/SX" 1.02 

for an 0(3) theory with monopoles in three dimensions and:- 

w*)> = ^ 

for an 0(2) theory with strings in three dimensions. The first observation is that the 
dependence of the density on time t is only through the correlation length ^{t). As the 
domains of coherent field form and expand, the inter defect distance grows accordingly. This 
we would interpret as the domains carrying the defects along with them on their boundaries. 
Secondly, there is roughly one defect per coherence size, a long held belief for whatever 
mechanism. However, in this case the density is exactly calculable. 

It is also possible to use Halperin's results to calculate defect-defect correlation functions. 
For example, the monopole-monopole correlation function on scales larger than a coherence 
length is found to be:- 

(Pnet(x)pnet(0))t = {p{t))tS{x) + ^(x)*, 

where g'(x)t is a measure of teh screening of a monopole at the origin. Explicit calculation 
yields:- 

3v^exp(-3a;Ve) sin3(2y2a;/0 



VI. CONCLUSIONS 



Under the conditions of a symmetry breaking phase transition, from 0{N) to 0{N — 1), 
which proceeds by a rapid quench, we have derived expressions for the evolution of the defect 
density during the early part of the fall from the false vacuum to the true vacuum. Our 
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results confirm that the defects are, indeed, frozen into the field when it first goes out of 
thermal equilibrium, at least for strongly coupled theories. Further, there is approximately 
one defect per correlation volume for the time during which the approximation is valid. The 
t^/^ time-dependence of the correlation-length ^(t) that we have seen above is specific to a 
non-conserved order parameter in a theory with a double time derivative, but we expect the 
qualitative features to be similar for all defect production by quenched symmetry breaking 
phase transitions. 

Thus, generically, we expect the defect density to follow the correlation length during 
the early part of the fall from the hill. Further, we expect the correlation length to grow 
more during this fall for weakly coupled theories, since the regime before which the fields 
feel the upturn in potential and slow down is longer. For very weakly coupled theories the 
correlation length can grow significantly in the time interval t = 0{ln{l/X) available. In 
almost all physically realistic scenarios except inflation, however, the coupling is not small 
enough for this growth to be significant. For more strongly coupled theories we know less 
Whatever, it is quite reasonable to say that the defect density is fixed at the time when the 
scalar field first goes out of equilibrium. 
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although see the self-consistent calulations of Boyanovsky et al. 



15 



REFERENCES 



[1] T.W.B. Kibble, J. Phys. A9, 1387 (1976). 

[2] T. Vachaspati and A. Vilenkin, Phys. Rev. D30, 2036 (1984). 

[3] V.I. Ginzburg, Fiz. Tverd. Tela 2, 2031 (1960); Sov. Phys. Solid State 2, 1826 (1961). 

[4] M. Hindmarsh and R.J. Rivers, Nud. Phys. 417, 506 (1994). 

[5] W.H. Zurek, Nature 317, 505 (1985). 

[6] W.H. Zurek, Acta Physica Polonica, B24, 1301 (1993). 

[7] B.I. Halperin, Les Houches, Session XXXV 1980 NATO ASI, editors Bahan, Kleman 
and Poirier. 

[8] K.T. Mahanthappa and P.M. Bakshi, J.M. Phys. 4, 1; ibid, 12 (1963) 
[9] J. Schwinger, J. Math. Phys. 2, 407 (1961). 
[10] L.V. Keldysh, Sov. Phys. JETP 20, 1018 (1965). 

[11] D. Boyanovsky, Da-Shin Lee and Anupam Singh, Phys. Rev. D48, 800 (1993). 
[12] D. Boyanovsky, H.J. de Vega, R. Holman. D-S. Lee and A. Singh, preprint DOR- 
ER/40682-77 



16 



